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SUMMARY 


^ The  optimization  problem  as  formulated  in  the  METRIC  model -f-^>  takes 
the  form  of  minimizing  the  expected  number  of  total  system  backorders 
in  any  two-echelon  inventory  system  subject  to  a budget  constraint. 

To  solve  this  problem,  one  needs  to  find  the  optimal  Lagrangian  multi- 
plier^associated  with  the  given  budget  constraint. 

For  any  large  scale  inventory  system,  this  task  is  computationally 
not  trivial.  Fox  and  Landi  proposed  one  method  which  was  a significant 
improvement  over  the  original  METRIC  algorithm.'")  In  this  report  we  first 
develop  a method  for  estimating  the  value  of  the  optimal  Lagrangian 
multiplier  used  in  the  Fox-Landi  algorithm,  present* alternative  ways 
for  determining  stock  levels,  and  sompar-a-,. these  proposed  approaches  with 

aSic  -c  i . 

the  Fox-Landi  algorithm^  using  two  hypothetical  inventory  systems — one 
involving  3 bases  and  75  items;  the  other  has  5 bases  and  125  items.  The 
comparison  shows  that  the  computational  time  can  be  reduced  by  nearly 
50  percent .— 

>Another  factor  that  contributes  to  the  higher  requirement  for  com- 
putational time  in  obtaining  the  solution  to  two-echelon  inventory  sys- 
tems is  that  it  has  to  optimally  allocate  stock  to  the  depot  as  well  as 
to  bases  for  a given  total  system  stock  level.  This  essentially  requires 
the  evaluation  of  every  possible  combination  o£\depot  and  base  stock 
levels — a time-consuming  process  for  many  practic^^inventory  problems 
with  a sizable  system  stock  level.  This  report  also  suggests  a simple 
approximation  method  for  estimating  the  optimal  depot  stock  level. 

When  this  method  was  applied  to  the  same  two  hypothetical  inventory 
systems  indicated  above,  it  was  found  that  the  estimate  of  optimal  depot 
stock  is  quite  close  to  the  optimal  value  in  all  cases.  Furthermore,  the 
increase  in  expected  system  backorders  using  the  estimated  depot  stock 
levels  rather  than  the  optimal  levels  is  generally  small. 


The  economic  interpretation  of ^a  Lagrangian  multiplier  is 
t ion  in  backorder  or  shortage  that  can  come  about  because  of  ar 
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Further  research  will  be  required  to  estimate  more  precisely  the 
reduction  in  computational  time  if  these  approximation  methods  are 
incorporated  in  the  requirements  computation  system  for  recoverable 
spares  (D041)  by  the  Air  Force.  For  instance,  we  would  have  to  apply 
the  proposed  methods  to  stratified  samples  of  recoverable  items  from 
the  D041  system  and  extrapolate  the  results  on  an  Air  Force-wide  basis. 
Nevertheless,  even  on  qualitative  grounds,  the  proposed  methods  are  so 
simple  and  reasonably  accurate,  our  conclusion  is  that  its  implementa- 
tion will  be  beneficial. 
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I.  introduction 


Almost  a decade  ago  Sherbrooke  formulated  the  well-known  METRIC 
model  for  determining  optimal  stock  levels  for  recoverable  items — 
items  subject  to  repair  when  they  fail — in  a two-echelon  setting  [3]. 
Briefly,  the  two-echelon  system  he  studied  consists  of  several  loca- 
tions, called  bases,  at  which  primary  demands  occur;  these  bases  are 
in  turn  resupplied  as  necessary  by  a central  repair  and  inventory 
stocking  point  called  a depot.  When  a failure  occurs  at  a base,  a 
demand  is  placed  on  base  supply  for  a corresponding  replacement  part. 
The  failed  part  is  then  either  repaired  at  that  base  or  is  sent  to  the 
depot  for  repair  depending  on  the  nature  of  the  failure.  Resupply  of 
base  supply  comes  from  the  base  maintenance  organization  if  repair  is 
accomplished  at  the  base;  otherwise,  resupply  comes  from  the  depot. 

In  either  case,  the  organization  resupplying  the  base  supply  activity 
does  so  by  exchanging  a serviceable  part  for  the  failed  part.  Thus  the 
inventory  policy  for  placing  orders  on  the  base's  maintenance  organiza- 
tion or  the  depot  is  an  (s  - l,s)  policy. 

Sherbrooke  presented  a model  for  determining  both  depot  and  base 
stock  levels  for  all  items  for  this  system.  In  particular,  the  problem 
he  formulated  was  to  minimize  the  average  total  number  of  base  back- 
orders existing  at  an  arbitrary  point  in  time  subject  to  a constraint 
on  system  investment,  that  is. 


m n 

min  £ l l (x 

1=1  i=l  x>s . . 

ij 


Sij )? (x I ^ijT 


ij 


(si0» 


(PI) 


m n 

subject  to  l £ c s < C , 
j=0  i=l  1J 
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= the  depot  stock  level  for  item  i, 

X = the  expected  daily  demand  rate  for  item  i at  base  j , 
c^  = the  unit  cost  for  item  i, 

C = the  budget  constraint, 

Tij  ^Si0^  = t*le  avera§e  resupply  time  for  base  j for  item  i 

given  the  depot  stock  level  for  item  i is  s^q,  and 
p(x|y)  = the  probability  that  x units  are  in  the  resupply 
system  given  that  the  expected  number  of  units  in 
the  resupply  system  is  y. 


Furthermore,  Sherbrooke  shows  that  T,.(s.„)  can  be  expressed  as 

ij  i0  K 


W * Vu + (1  - V(Bij  - 5<si0>  • V • 


A„  = the  average  base  repair  time  for  item  i at  base  j , 
r_  = the  proportion  of  demands  requiring  base  repair 
for  item  i at  base  j , 

B„  = the  average  order-and-ship  time  at  base  j for 
item  i, 

= the  average  depot  repair  cycle  time  for  item  i, 
6(s.q)  ’ = (1/X1)  l (x  - s1Q)p(x|XiDi) , the  expected 

X>SiO 

delay  per  depot  demand  for  item  i,  and 


L = £ (1  - t*ie  expected  daily  depot  demand 


for  item  i. 


In  the  remainder  of  the  report,  i will  refer  to  an  item  and  j will 
refer  to  a base  (j  = 0 represents  the  depot).  Thus  i and  j will  al- 
ways be  elements  of  the  sets  { 1 , n}  and  { 0 , m} , respectively. 

Additionally,  an  integer  k appearing  in  the  text  to  the  right  of  the 
statement  of  a problem  or  equation  will  designate  for  future  reference 
that  problem  or  equation.  For  complete  description  of  problem  back- 
ground and  formulation,  see  Ref.  3. 
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Subsequently  Fox  and  Landi  suggested  a Lagrangian  approach  for 
solving  problem  PI  [2].  One  obstacle  to  the  implementation  of  METRIC 
using  the  Fox-Landi  algorithm  is  the  requirement  of  estimating  an  ap- 
propriate value  for  the  Lagrangian  multiplier.  Another  important  and 
related  problem  is  the  lengthy  computer  run  time  required  to  obtain 
an  optimal  solution  to  problem  PI  when  using  their  algorithm.  A large 
portion  of  this  computational  effort  is  related  to  searching  for  the 
optimal  depot  stock  level.  This  search  is  particularly  time-consuming 
for  items  having  a high  average  number  of  units  in  the  depot  repair 
cycle  since  the  amount  of  computation  required  by  their  algorithm  is 
roughly  proportional  to  the  number  of  depot  stock  levels  explicitly 
examined . 

In  this  report  we  first  develop  an  approach  for  obtaining  an  esti- 
mate of  the  optimal  Lagrange  multiplier  value  required  in  the  Fox-Landi 
algorithm,  present  two  new  methods  for  determining  stock  levels,  and 
compare  these  methods  with  the  Fox-Landi  method  and  other  techniques. 
The  proposed  approach  eliminates  the  particularly  time-consuming  por- 
tion of  the  Fox-Landi  algorithm  devoted  to  searching  for  the  best 
Lagrange  multiplier  value  and  significantly  reduces  computation  time 
for  determining  stock  levels  without  degrading  the  quality  of  the 
solution. 

We  then  present  a method  for  estimating  the  optimal  depot  stock 
level.  Limited  computational  experience  indicates  that  this  method  is 
easy  to  implement,  provides  a very  good  estimate  of  the  optimal  depot 
stock  level,  and  is  particularly  useful  for  items  having  a high  average 
number  of  units  in  the  depot  repair  cycle.  For  these  items  it  is  pos- 
sible to  reduce  computation  time  required  by  the  Fox-Landi  algorithm 
by  as  much  as  90  percent. 


f 
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II.  THE  APPROXIMATION  PROBLEM 


In  this  section  we  first  construct  a problem  that  is  a continuous 
approximation  to  problem  PI.  We  then  state  and  prove  two  theorems  that 
are  the  basis  for  an  algorithm  that  can  be  used  to  solve  this  approxi- 
mating problem. 

Recall  that  the  total  average  base  backorders  existing  at  any 
point  in  time  for  item  i can  be  expressed  as 


m 

l l 

j=l  x>s . . 

ij 


(x 


Sij)P(xlXijTij(siO)) 


Two  useful  probability  distributions  for  describing  the  demand  process 

are  the  Poisson  and  negative  binomial  distributions.  As  shown  in 

Ref.  1,  this  implies  that  if  demand  has  a Poisson  or  negative  binomial 

distribution,  then  for  a given  value  of  X . .T . . (s  ._) , p (x I X . .T . . (s . _) ) , 

i]  i]  lO  K 1 ij  ij  lO 

the  probability  distribution  representing  the  number  of  units  in  re- 
supply of  item  i at  base  j at  any  point  in  time,  is  a Poisson  or  nega- 
tive binomial  distribution,  respectively. 

Experimental  data  gathered  during  the  conduct  of  this  study  indi- 
cate that  when  p (x | X. .T . . (s .„) ) is  either  a Poisson  or  negative  binomial 
1 lj  lj  lO  “ 

distribution,  the  above  total  expected  backorder  expression  can  be 
closely  approximated  by  an  exponential  function.  That  an  exponential 
function  accurately  approximates  this  expression  should  not  be  entirely 
unexpected.  First,  for  budgets  of  practical  interest,  the  item  stock 
levels,  s^j>  are  normally  much  larger  than  the  average  demand  during 
the  resupply  time.  In  fact,  the  probability  of  running  out  of  stock 
during  the  resupply  time  is  often  much  less  than  .15  in  real  applica- 
tions. Thus  the  only  probabilities  entering  the  backorder  calculation 
are  the  tail  probabilities  of  the  distribution.  In  the  tails,  both  the 
Poisson  and  negative  binomial  distributions  behave  almost  like  the 
geometric  distribution;  that  is,  each  succeeding  probability  is  roughly 
a constant  proportion  of  its  predecessor.  Consequently,  when  s^  is 


large  relative  to  ^si0^ 4 t*le  exPecte<^  number  of  backorders 

existing  at  any  time  at  location  j for  item  i is  approximately  a geo- 
metric function  of  s„  . Therefore,  an  exponential  function  is  a use- 
ful continuous  approximation  to  this  relationship  between  expected 
backorders  at  a location  and  the  item's  stock  level  at  that  location. 

Furthermore,  total  expected  base  backorders  exhibit  this  same 
behavior.  If  demand  has  either  a Poisson  or  negative  binomial  dis- 
tribution (or,  for  that  matter,  any  compound  Poisson  distribution), 
then  the  total  number  of  units  of  an  item  in  resupply  across  all  bases 
has  a Poisson  or  negative  binomial  distribution,  respectively,  given 
we  assume  independence  of  demand  and  common  var iance-to-mean  ratio 
among  bases.  Since  in  most  practical  situations  total  system  stock 
substantially  exceeds  the  total  expected  number  of  units  in  resupply, 
the  tail  of  the  distribution  describing  the  total  number  of  units  in 
resupply  is  the  only  portion  of  the  distribution  of  importance.  As  an 
approximation,  this  distribution  can  be  used  to  determine  the  nature  of 
the  relationship  between  total  expected  base  backorders  and  total  system 
stock.  For  the  reasons  discussed  previously,  an  exponential  function 
should  also  adequately  represent  this  relationship. 

Thus  we  will  approximate 


I l <* 

jBl  *>sij 


s . . )p (x  A . . T . . (s 
XJ  ij 


iO 


)) 


with  the  exponential  function 


-b.N. 

B.  (N.)  = a.e  11 
11  l 

In  this  approximation,  N.  represents  total  system  stock.  In  practice, 
the  parameters  > 0 and  b.  > 0 are  estimated  using  regression  analy- 
sis. The  data  used  in  the  regression  analysis  are  the  backorder  data 
obtained  from  the  solution  to  the  problem 
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min  I l (x  - s )p(x|A  T (s  )) 

i=i  x>s . . J j j 

J 1 -1 


subject  to 


s . . = N . , and 
ij  i 


s. . =0,  1,  . . . , N.  , 
xj  1 


for  several  appropriate  values  of  N^. 

We  now  formulate  a continuous  approximation  to  problem  PI  in 
which  the  exponential  representation  of  total  system  backorders  for 
an  item  is  used.  In  this  approximation  problem,  the  decision  vari- 
ables are  the  total  system  stock,  N^,  rather  than  the  stock  levels 
for  each  location,  s„.  As  we  shall  see,  the  main  reason  for  studying 
this  approximation  problem  is  that  it  is  a vehicle  for  obtaining  an 
estimate  of  the  optimal  Lagrangian  multiplier  value  used  in  the  Fox- 
Landi  algorithm.  The  approximation  problem  is  formulated  as  problem  P2: 


min  V B (N . ) 
i=l  1 1 


subject  to  £ c.N.  < C , 
i=l 


where  > 0 


Note  that  is  a continuous  variable  in  this  approximation.  The 
optimality  conditions  (Kuhn-Tucker  conditions)  for  this  problem  are 
as  follows: 

Find  0^  > 0 such  that 


<a)  dMj  + 6lCi  i ° ■ 


(b)  l c N < C , 
i=l 


N.  > 0 , 


(c)  9i  i/A'T0’ 


(d>  Ni  l dir  + 0lCi)  = 0 • 


A relaxed  version  of  problem  P2  in  which  the  non-negativity  con- 
straint on  the  item  stock  level  is  removed  is  problem  P3: 


lin  I B . (N . ) 
i=l  1 


subject  to  ][  ci^i  = ^ 
i=l 


The  optimality  conditions  for  this  problem  are: 
Find  ©2  ^ 0 such  that 


(a)  diT  + 92Ci  = ° • 


(b)  l c N < C , 
i=l 


(C)  02  ( CiNi  - C)  = ° * 


(d)  Ni  ( diTt  + 02Ci)  = ° • 


UK's  w?  ' 
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We  now  explore  the  relationship  between  problems  P2  and  P3  in 
detail . 

* 1 

Suppose  we  obtained  a solution  to  problem  P3.  Let  N.  represent 

2 1 

the  optimal  solution  to  problem  P2,  and  N.  represent  the  optimal  solu- 

2 1 1 2 

tion  to  problem  P3.  If  N.  >0  for  all  i,  then  N.  = N.  and  the  obiec- 

i = it  J 

tive  function  values  are  equal. 

2 

Suppose,  however,  that  < 0 for  at  least  one  value  of  i.  Let 

= max(0,  N^) 


and 


n 

= I 

i=l 


c ,N . 
i l 


— 2 „ — 2 — 

Since  N.  > N.  for  all  i and  N.  > N.  for  at  least  one  i,  C > C. 

1=1  ii  ’ 

Suppose  problem  P2  is  modified  slightly  so  that  the  right-hand 
side  value  C is  replaced  by  C.  This  modified  problem  is  problem  P4: 


n 

min  y B . (N  . ) 
i=l 


n 

subject  to  £ c.N.  C , (P4) 

i=l  1 1 

where  > 0 . 

The  optimality  conditions  for  this  problem  are  the  same  as  those  given 
for  problem  P2  after  substituting  C for  C.  Also,  let  0 represent  the 
optimal  value  of  the  Lagrangian  multiplier  for  problem  P4 . 

In  solving  problem  P3,  we  will  obtain  a value  for  0„.  We  now 

— _ 2 ^ 
show  that  0 = 0£,  and  that  = max(0,  N^)  is  an  optimal  solution  to 

* 

Section  III  develops  the  method  for  determining  the  solution  to 
problem  P3. 


w 
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problem  P4  by  demonstrating  that  these  values  satisfy  the  Kuhn-Tucker 
conditions  corresponding  to  problem  P4. 

By  construction. 


ri  _ _ _ _ / n 

y c.N.  = C , N.  > 0 , and  6 ) c.N.  - C I = 0 

i=l  11  1 \i=l  1 1 ' 


If  0 = Q^,  9^0  since  0^  > 0.  Suppose  N^  = N'2,  that  is,  N2  > 0. 
Then 


dB. 

l 

dN. 

l 


dB. 

l 

dN. 


N.=N. 

l l 


2 

N.=NT 

1 l 


and 


0 = 


dB. 
i 

dN. 

l 


dB. 

„ + 92Ci  = dN? 
N.=N2 

l l 


+ 0c.  . 

l 


N=N . 

l 


By  assumption  there  exists  at  least  one  value  of  i for  which 

— 2 — 2 
N.  > N. ; that  is,  N.  = 0 while  N.  < 0.  Since 
ii  i i 


dB. 

dN. 


dB. 

l 

dN. 


N . =0 

l 


N <0 

l 


due  to  the  exponential  form  of  B^(N^),  and 


dB. 

l 

dN. 


2 

N =N; 
1 1 


+ 02ci  = ® > 


we  know  that 
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dN. 


+ 6c  > 0 . 


N1=0 


Consequently,  the  optimal  solution  to  problem  P4  is  N.  = N.  = 

2 1 
max{0,  N^}.  Furthermore,  the  optimality  conditions  are  satisfied 

when  0 is  equal  to  02- 

Theorem  1 . 0^  ©2- 

Proof;  The  optimal  objective  function  value  for  problem  P2  is  a 
convex,  differentiable,  strictly  decreasing  function  of  the  available 
budget,  C.  Since  the  slope  of  this  function  is  equal  to  the  negative 
of  the  Lagrangian  multiplier  value,  0^  £ 0 since  C £ C.  But  = 0, 
so  0^  ©2* 

Corollary.  0^  > ©2  when  C > C. 


Next  we  compare  N^  with  N^. 


If  C = C,  then  N^  = N.  for  all  i. 
l i 


Now  let  us  suppose  C > C so  that  0^  > ©2 


0.  Let  us  examine  the  two 


cases  N^  > 0 and  N^ 


First,  assume  N^  > 0 


0 separately. 
Then 


d^ 

dN. 


+ 0c^  = 0 


Ni=Ni 


Furthermore,  if  N^  > 0,  then 


d^ 

dN, 


+ 0lCi  = 0 


elci  > 9ci 


dB 

dN. 


VNi 


Since 
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In  this  section  we  established  several  important  relationships 
among  problems  P2,  P3,  and  P4.  In  Sec.  Ill  we  develop  a simple  algo- 
rithm for  solving  problem  P2  based  on  these  relationships  and  show 
how  to  find  the  solution  to  problem  P3.  As  we  have  just  demonstrated, 
once  we  have  the  solution  to  problem  P3,  we  also  have  the  solution  to 
problem  P4.  From  Theorem  2,  we  then  have  an  upper  bound  on  the  values 
of  the  N^.  In  particular,  if  = 0,  then  = 0.  Combining  this 
observation  with  the  implications  of  Theorem  1 and  its  corollary  pro- 
vides the  basis  for  the  proposed  algorithm  for  solving  problem  P2. 
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III.  COMPUTING  OPTIMAL  SOLUTIONS  FOR  PROBLEMS  P2  AND  P3 


We  begin  this  section  by  developing  a method  for  determining  the 
optimal  solution  to  problem  P3.  Observe  that  the  optimal  solution 
must  satisfy  the  following  two  conditions: 


dB. 

sr + Vi  ■ 0 

1 


and 


I C N.  = C 
i-1  1 1 


The  second  condition  must  hold  since  each  B^(N^)  is  a strictly  de- 
creasing function  of  N^. . 

Since 


B.(N.)  = aie 


-biNi 


where  a^,  b^  > 0,  the  first  condition  states  that 


-biNi 


> 0 , 


or 


'Vi 


d 


i 


Letting 


where  g^  = - a and  f^  = Bb^.  Consequently,  IT  is  a linear"  function 

of  C.  If  the  budget  is  incremented  by  an  amount  AC,  then  N!,  the  new 
value  of  the  stock  level  for  item  1,  satisfies 


, AC 

N'  = N + — 
i if* 


The  optimal  solution  to  problem  P2  has  been  found  if  each  of  the 
found  using  Eq.  (E2)  is  non-negative.  If  there  exists  an  i for 

which  N.  < 0,  then  we  may  employ  the  following  algorithm  to  find  the 

* 1 
optimal  solution  to  problem  P2.  Let  I = {l,  ...,  n}  and  represent 

the  optimal  solution  to  problem  P2. 

Step  0.  Solve  Problem  P3  as  described  above,  thereby  obtaining 

an  initial  value  for  N , i e I. 

1 1 

Step  1.  Set  = 0 for  all  < 0 during  the  last  iteration  and 
delete  the  corresponding  i from  I.  Recompute  a and  B,  where 


a-  l 

iel 


r 

iel  I i 


Step  2.  Using  Eq.  (E2),  obtain  new  estimates  of  for  each 

iel.  If  ^ 0 for  all  iel,  then  the  optimal  solution  has  been 

found,  and  = N.  for  all  iel  and  N'f  = 0 for  all  i = 1,  ...,  n for 

i i l 

which  i £ I.  If  there  exists  some  i for  which  N.  <0,  return  to 
step  1. 

It  is  clear  that  our  solution  satisfies  all  the  optimality  condi- 
tions for  problem  P2  except  possibly  condition  (a)  for  i £ I.  However, 
at  an  earlier  iteration  (when  i was  deleted  from  I)  we  had 
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= 0 


where  "8  and  {^(<0)  are  the  earlier  values  of  and  N.,  respectively. 
Since  dBj/dN^  is  clearly  increasing  in  N^,  and  02  increases  at  each 
iteration  (Theorem  1 and  its  corollary),  condition  (a)  must  hold.  Con 
vergence  is  guaranteed  since  n is  finite. 
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IV.  A COMPARISON  OF  ALTERNATIVE  SOLUTION  PROCEDURES 
FOR  SOLVING  PROBLEM  PI 


In  this  section  we  review  three  algorithms  for  solving  problem  PI 
and  compare  them  to  two  algorithms  designed  to  obtain  a solution  for 
problem  PI  based  on  the  solution  to  the  approximating  problem,  problem 
P2. 

THE  SHERBROOKE  PROCEDURE 

The  first  algorithm,  a procedure  originally  proposed  by  Sherbrooke 
[3],  is  a marginal  analysis  algorithm  consisting  of  two  phases.  In  the 
first  phase,  each  item  is  examined  independently.  The  optimization 
problem  solved  for  item  i in  the  first  phase  has  the  form: 

m 

min  l l (x  - s )p(x|X  T (s  )) 

J=1  x>s..  J J 


subject  to 


3=0 


13 


= N. 


where 


1,  • 


(P5) 


and  N.  is  the  total  system  stock  available  for  distribution  among  the 
depot  and  bases.  Let  Z (N^)  represent  the  optimal  objective  function 
value  for  problem  P5  given  units  are  available  for  distribution. 
Problem  P5  is  solved  by  obtaining  the  solution  to  the  N + 1 problems 


^i(Ni’  Si0)  = min  j_  (x  ” sij)p(x|AijTij(sio)) 


j=l  x>s 


13 


subject  to  £ s = - s^  , 


J-l 


(P6) 


where  s^  =0,  1,  ...» 
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and  is  fixed  for  =0,  1,  N^.  Problem  P6  can  be  solved 

via  marginal  analysis . Then 

Z.(N.)  = min  Z^N.,  s^)  , 
si0 

where  s^q  =0,  . . . , . 

The  second-phase  problem  is 

n 

min  T Z . (N . ) 
i=l  1 1 

n 

subject  to  £ ciNi  = ^ » 
i=l 

where  N.  =0,  1,  ...  . 

l 

Sherbrooke  [3]  suggests  that  a marginal  analysis  algorithm  be  used  to 
find  a solution  to  this  knapsack  problem.  Clearly  other  procedures 
could  be  employed  to  obtain  an  optimal  solution.  In  any  case,  this 
approach  requires  a substantial  amount  of  storage  to  save  all  the 
Z (N.)  values.  For  moderately  sized  problems — several  thousand  items — 
a storage  requirement  of  10  or  more  words  may  be  needed  to  save  these 
values.  Furthermore,  the  computation  time  required  to  obtain  these 
^i<^i>  va^ues  ^or  such  problems  is  very  large. 

THE  FOX  AND  LAND I PROCEDURE 

Subsequently  Fox  and  Landi  [2]  proposed  a Lagrangian  algorithm 
for  solving  problem  PI.  In  particular,  they  formulated  the  relaxed 
version  of  problem  PI  as  problem  P7: 
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m n 

min  l l I (x 

1=1  i=l  x>s , . 

ij 


where 


s,  . 
ij 


0,  1, 


Sij)P(xlAijTij(sio))  + 0 


I I 

j=0  1=1 


c.s.  . 
i ij 


(P7) 


and  0 is  the  Lagrangian  multiplier.  Since  problem  P7  is  separable  by 
item,  its  optimal  solution  can  be  found  by  solving  the  n individual 
item  problems 


"1"  A is  <x ' V’^Wio”  + 6 J0  cisu 


subject  to  s^  =0,  1,  ... 


This  problem,  like  problem  P6  in  Sherbrooke's  two-phase  method,  is 
solved  using  a partitioning  procedure,  that  is,  it  is  reformulated  as 


min  < 0c.s._  + 7 min 

^ , i i0 


siQ  0,1,... 


j=l  s1j=0 , 1 , . 


J (x  ' sij)p(xlXijTij(sio))  + 0cisij:  sio  fixedi 
X Sij  ' 


\ , 


(P8) 


or  equivalently  as 


min  Z(siQ;  0) 

where  s^  = 0,  1 , ...  , 


(P9) 


and 
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Z(s.0;  0)  = ec.s10  + (x  - s. . )P (x 1 A. . T . . (s.0) ) 


1=1  s . . I x>s . . 
1J  v 11 


+0C.S..:  s..  =0,  1,  ...  ; s._  fixed 
l ij  ij  lO 


I. 


To  determine  Z(s^q;  0),  solve  the  m base  problems 


min  l (x  - s..)p(x|A..T..(s  ))  + 0c. s.. 

s..  x>s  1J  10  1 ^ 

iJ  ij 


The  optimal  s is  the  smallest  non-negative  integer  for  which 


l P(XU  T..(Si0))  < 0ci  • 


x>s . . 

ij 


Problem  P8  is  solved  for  each  item  for  a given  value  of  0.  This 
yields  a total  investment  cost  corresponding  to  0.  In  the  Fox-Landi 
approach,  the  "optimal"  value  of  0 is  selected  from  a grid  of  M equally 
spaced  values. 


% > 01  > •••  > »„  > 0 • 

The  optimal  value  of  0 is  the  0V,  K £ {0,  ...,  M},  whose  corresponding 

K. 

total  investment  cost  is  closest  to  C. 

Fox  and  Land!  suggest  that  their  method  is  a single-pass  method; 
that  is,  only  one  pass  through  the  item  data  base  is  necessary  to  ob- 
tain the  optimal  solution.  The  storage  requirement  to  effect  this 
one-pass  approach  is  potentially  enormous.  For  a moderately  sized 
problem  having  3000  items,  20  bases,  and  M = 63,  almost  4 million  item 
stock  levels  must  be  saved,  plus  possibly  millions  of  additional  item 
data  elements  reflecting  fill  rates,  probability  of  no  stockout  at  an 
arbitrary  point  in  time,  expected  base  backorders,  etc.  Furthermore, 
because  there  may  be  no  simple  method  for  estimating  suitable  bounds 


1 


A 
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on  the  values  of  the  multipliers,  much  larger  values  of  M may  be  re- 
quired to  ensure  adequate  approximation  of  the  budget. 

It  has  been  the  author's  experience  that  Air  Force  personnel  have 
difficulty  estimating  a reasonable  range  for  0 for  large  problems. 

This  is  not  surprising  because  the  data  used  in  the  model  frequently 
change  in  real  situations,  thereby  causing  the  optimal  value  of  the 
multiplier  to  change.  Furthermore,  changing  the  multiplier's  magnitude 
by  10  ^ or  less  often  causes  the  corresponding  total  cost  to  change  by 
many  millions  of  dollars.  Consequently,  2^  values  of  0 have  been  used 
in  some  Air  Force  applications  to  make  the  system  "foolproof."  In  these 
cases  60  million  or  more  item  stock  levels  would  be  needed  to  be  ex- 
plicitly stored — plus  a considerable  amount  of  other  item  and  base 
data — to  make  the  Fox-Landi  algorithm  truly  a one-pass  method. 

On  the  other  hand,  if  their  method  is  altered  so  that  the  item 

data  are  passed  through  a second  time,  it  is  possible  to  eliminate 

virtually  all  the  requirement  for  secondary  storage,  In  the  first 

pass,  only  the  running  total  cost  corresponding  to  each  0,., 

K. 

K £ {0,  ...,  M},  is  saved.  At  the  end  of  this  phase  the  "optimal" 
multiplier  value,  0*,  is  established.  The  second  phase  of  the  algo- 
rithm requires  a second  pass  through  the  data  base.  In  the  second 
pass,  the  optimal  stock  levels  for  each  location  are  found  for  all 
items  by  resolving  problem  P8  with  0 = 0*. 

In  some  applications  the  Fox-Landi  one-pass  method  is  clearly 
infeasible;  that  is,  there  may  not  be  enough  peripheral  storage  ca- 
pacity to  save  all  the  data.  If  storage  capacity  is  available,  there 
is  a tradeoff  between  the  time  and  cost  required  to  store  and  access 
the  data  in  the  secondary  memory  using  the  one-pass  method,  and  the 
time  and  cost  to  recompute  the  stock  levels  using  the  second  method. 

For  realistic  Air  Force  problems,  the  two-pass  method  appears  to  be 
the  only  feasible  approach  given  current  hardware  constraints  if  M is 
large  enough  to  guarantee  that  a solution  can  be  found  that  closely 
approximates  the  target  budget. 
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THE  BISECTION  METHOD 

A third  way  to  solve  problem  PI  is  a slight  modification  for  the 
Fox-Landi  algorithm  called  the  bisection  method,  which  employs  a bi- 
section search  to  find  the  optimal  value  for  0.  This  procedure  requires 
initial  upper  and  lower  bounds  on  the  optimal  value  of  0.  Call  these 
0jj  and  0^,  respectively.  The  bisection  method  is  as  follows: 

1.  Set  0 = (Oy  + 0l)/2. 

2.  Solve  problem  P8  with  0=0  for  each  item. 

3.  If  the  total  cost  of  the  solution  obtained  in  step  2 exceeds 
C,  then  replace  0^  with  0;  otherwise,  replace  0 with  0‘ 

4.  If  a stopping  criterion  has  not  been  met  (such  as  a fixed 
number  of  iterations  or  an  error  tolerance),  return  to  step 
1;  otherwise,  stop. 

The  major  drawback  to  the  bisection  approach  is  that  a separate 
pass  through  the  item  data  base  is  required  at  each  iteration  of  the 
algorithm.  This  algorithm  performs  very  well  in  terms  of  convergence, 
and  we  have  found  that  it  almost  always  produces  solutions  that  are 
within  1/2  percent  of  the  target  budget  using  10  bisections. 

COMPARISON  OF  METHODS 

The  closeness  of  the  solutions  to  the  target  budget  generated  by 
either  the  Fox-Landi  method  or  the  bisection  algorithm  depends  on  how 
broad  a range  of  multiplier  values  must  be  searched  for  a fixed  value 
of  M or  a fixed  number  of  bisections.  It  should  be  pointed  out  that 
both  of  these  methods  only  yield  an  approximation  to  the  optimal  multi- 
plier value  (assuming  one  exists). 

Of  the  methods  discussed  thus  far,  it  has  been  the  experience  of 
the  author,  as  well  as  of  Fox  and  Landi  [2],  that  the  latter  two  algo- 
rithms dominate  Sherbrooke's  algorithm  in  run  times  by  an  order  of  mag- 
nitude or  more  on  real  problems  given  reasonable  estimates  of  upper  and 
lower  bounds  for  the  Lagrangian  multiplier.  Thus  in  the  comparisons 
we  will  report,  only  these  two  Lagrangian  methods  will  be  discussed. 
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APPROXIMATION  METHODS 

Earlier  we  described  an  approximation  method  for  estimating  the 
optimal  values  of  0 and  each  N^.  Several  options  are  open  for  imple- 
menting this  approximation  method.  One  way  to  implement  it  is  to  use 
a two-phase  approach.  Call  this  approach  the  First  Approximation 
Method.  The  values  of  a^  and  b^  are  computed  in  the  first  phase  of 
this  method  and  the  optimal  value  of  0 is  estimated  using  Eq . (El). 

In  the  second  phase,  we  solve  problem  P8  for  each  item,  using  the  esti- 
mate of  the  optimal  0.  This  approach  has  two  major  advantages  over  the 
Fox-Landi  method: 

1.  The  estimate  of  the  optimal  multiplier  can  be  obtained  without 
prespecifying  a range  of  values,  and  computation  time  to  ob- 
tain the  estimate  does  not  depend  on  the  uncertainty  of  the 
multiplier  value. 

2.  The  computation  time  to  find  an  estimate  of  the  optimal  multi- 
plier is  much  smaller. 

If  the  two-pass  version  of  the  Fox-Landi  algorithm  is  used , the  second 
phase  of  that  method  and  the  approximation  method  are  the  same.  The 
one-pass  version  of  the  Fox-Landi  algorithm  requires  considerably  more 
storage,  and  also  requires  more  computer  time  to  determine  the  optimal 
stock  levels  than  this  approximation  method  requires. 

This  approximation  approach  also  has  advantages  over  the  bisection 
method : 

1.  Only  two  passes  through  the  data  base  are  required  as  opposed 
to  seven  or  more  required  for  the  bisection  method  in  practice. 

2.  No  stock  levels  need  to  be  saved;  in  the  bisection  method  it 
is  necessary  to  save  all  stock  levels  and  other  data  for  three 
multiplier  values. 

Another  algorithm  can  be  employed  that  directly  uses  the  results 
of  the  approximation  problem,  that  is,  problem  P2.  Call  this  approach 
the  Second  Approximation  Method.  This  algorithm  is  of  interest  in 
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situations  in  which  we  only  want  to  compute  total  system  stock  for 
each  item  and  are  not  particularly  interested  in  computing  the  optimal 
distribution  of  the  assets.  Determining  the  optimal  allocation  of  a 
budget  among  items  is  of  primary  importance  when  purchasing  inventory 
or  making  budgetary  projections  for  spares  for  different  systems.  In 
these  cases,  distribution  decisions  are  usually  not  that  critical. 

This  Second  Approximation  algorithm  also  consists  of  two  phases. 

In  the  first  phase  we  estimate  the  values  of  the  a.  and  b.  parameters, 

1 1 

and  in  the  second  phase  we  determine  total  system  stock  for  each  item 
using  the  algorithm  described  in  Sec.  Ill  and  rounding  N.  to  the 
nearest  integer.  The  algorithm  requires  one  pass  through  the  item 
data  base  and  one  pass  through  an  item  file  consisting  of  a^,  h.,  and 
c^.  The  major  advantage  of  this  approach  is  that  it  eliminates  the 
stock  allocation  phase  of  both  the  Fox-Landi  method  and  the  First 
Approximation  algorithm. 
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V.  A COMPUTATIONAL  COMPARISON  OF  VARIOUS  ALGORITHMS 


The  Fox-Landi  algorithm,  bisection  algorithm,  and  the  two  approxi- 
mation methods  have  been  coded  and  tested  on  several  sample  sets  of 
data  for  the  Air  Force's  new  F-15  fighter.  Since  all  of  the  tests 
yielded  the  same  general  results,  we  will  discuss  only  two  of  them 
in  detail.  The  first  test  consisted  of  a 75-item  sample  and  had  3 
operating  bases.  The  flying  programs  were  very  different  at  each  base. 

In  the  second  test,  125  items  were  included  in  the  sample  with  demands 
occurring  at  5 bases.  In  the  second  test,  only  the  Fox-Landi  and  the 
two  approximation  methods  were  compared.  The  run  times  stated  for  both 
approximation  algorithms  include  the  time  required  to  estimate  the 

values  of  a.  and  b..  In  all  Fox-Landi  calculations,  a maximum  of  128 
l l 

multiplier  values  were  examined;  ten  bisections  were  used  in  all  appli- 
cations of  the  bisection  method.  Furthermore,  in  both  test  cases  all 
stock  levels  for  all  relevant  multiplier  values  were  stored  in  main 
memory.  Thus,  although  the  reported  computation  times,  which  include 
compile  times,  are  roughly  equal  for  all  the  algorithms,  they  are 
biased  in  favor  of  the  Fox-Landi  method  because  this  type  of  storage 
would  be  impossible  for  larger  problems.  In  addition,  the  range  of 
multiplier  values  considered  in  the  test  of  the  Fox-Landi  and  bisection 
methods  was  selected  after  estimating  the  optimal  multiplier  value 
using  the  First  Approximation  Method.  Thus  the  test  results  are  biased 
in  favor  of  them,  since  the  range  of  multiplier  values  was  much  smaller 
than  would  normally  be  the  case. 

The  data  displayed  in  Tables  1 and  2 indicate  how  well  each  approach 
approximates  a given  target  budget  for  the  two  test  data  sets.  Without 
a doubt  the  bisection  method  produced  solutions  that  best  matched  the 
target  budgets,  followed  in  order  by  the  Second  Approximation  Method, 
the  Fox-Landi  method,  and  the  First  Approximation  Method.  As  mentioned 
before,  the  results  are  biased  in  favor  of  both  the  Fox-Landi  and  bi- 
section methods  due  to  the  initialization  of  the  range  of  multiplier 
values.  From  a practical  viewpoint,  all  approaches  worked  acceptably 
well  in  meeting  the  target  budgets.  Furthermore,  the  stock  levels 


Table  1 


7 5- ITEM,  3-BASE  TEST  CASE 


Total  Cost  ($  millions) 

Target 

Budget 

Bisection 

Fox-L  .ndi 

First 

Approximation 

Second 

Approximation 

3.68 

3.67 

3.68 

3.63 

3.63 

3.97 

3.99 

3.92 

3.82 

4.03 

4.27 

4.27 

4.27 

4.30 

4.18 

4.57 

4.57 

4.57 

4.62 

4.61 

4.87 

4.87 

4.85 

4.87 

4.78 

5.16 

5.16 

5.18 

5.09 

5.17 

5.46 

5.46 

5.42 

5.38 

5.49 

5.76 

5.76 

5.76 

5.75 

5.79 

6.05 

6.06 

6.05 

6.06 

6.08 

6.35 

6.34 

6.28 

6.28 

6.33 

6.65 

6.65 

6.63 

6.63 

6.73 

6.94 

6.89 

6.80 

6.  87 

6.92 

7.24 

7.24 

7.19 

7.27 

7.24 

7.54 

7.54 

7.57 

7.68 

7.51 

7.83 

7.84 

7.77 

7.80 

7.83 

8.13 

8.14 

8.24 

8.20 

8.05 

8.43 

8.42 

8.50 

8.42 

8.42 

8.73 

8.73 

8.50 

8.74 

8.77 

9.02 

9.02 

9.04 

9.11 

9.00 

Execut ion 
time 
(sec) 

92.57 

19.57 

11.59 

4.57 

generated  by  the  various  approaches  were  virtually  the  same  for  similar 
budgets.  Consequently,  total  system  expected  backorders,  for  all  prac- 
tical purposes,  are  indistinguishable;  that  is,  the  backorder  versus 
investment  curves  virtually  coincide  among  these  various  approaches. 
Exact  comparison  of  computed  stock  levels  and  expected  backorders  can- 
not be  made  among  the  competing  methods  since  the  allocation  of  the 
available  budget  in  each  case  depends  on  the  way  each  algorithm  esti- 
mates the  Lagrangian  multiplier. 

The  area  in  which  the  methods  clearly  differ  is  in  computation 
time.  The  approximation  methods  require  substantially  less  time  than 
either  the  Fox-Landi  method  or  the  more  time-consuming  bisection  method. 
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Table  2 

125-ITEM,  5-BASE  TEST  CASE 


Total  Cost  ($  millions) 


Target 

Budget 

Fox-Landi 

First 

Approximat ion 

Second 

Approximation 

26.4 

26.7 

24.8 

26.6 

27.6 

27.6 

26.2 

27.9 

28.7 

28.7 

27.6 

28.9 

29.8 

30.0 

29.5 

29.8 

31.0 

31.2 

30.7 

30.8 

32.1 

32.1 

32.0 

32.2 

33.2 

33.3 

33.1 

33.1 

34.4 

34.4 

34.3 

34.2 

35.4 

35.5 

35.9 

35.7 

36.6 

36.8 

37.0 

36.7 

37.8 

38.0 

38.1 

37.7 

38.9 

38.6 

39.3 

39.2 

40.0 

39.9 

40.6 

40.0 

41.2 

41.1 

42.1 

41.3 

42.3 

42.5 

43.9 

42.4 

43.4 

43.3 

44.7 

43.7 

44.6 

44.5 

45.6 

44.2 

45.7 

46.3 

46.1 

45.9 

46.8 

47.2 

47.3 

46.7 

Execution 

time 

(sec) 

36.98 

16.28 

4.74 

NOTE:  All  programs  are  run  on  an  IBM  370/168. 


Other  experimentation  has  shown  that  the  percentage  difference  in  com- 
putation times  tends  to  be  even  more  substantial  as  the  number  of  items 
considered  increases. 

Thus  the  approximation  methods  produce  answers  that  are  as  good  as 
those  produced  by  the  Fox-Landi  method  and  the  bisection  method,  only 
much  more  quickly  than  those  methods.  The  bisection  method  does,  how- 
ever, match  target  budgets  slightly  better  than  the  approximation 
methods.  However,  the  approximation  algorithms  are  virtually  foolproof, 
which  is  perhaps  their  greatest  advantage.  The  user  does  not  have  to 
specify  the  range  of  multiplier  values  or  the  number  of  bisections  in 


advance.  This  eliminates  one  problem  associated  with  implementing 
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VI.  ESTIMATION  OF  THE  OPTIMAL  DEPOT  STOCK  LEVEL 


We  have  described  Sherbrooke's  algorithm  and  several  Lagrangian 
type  methods  for  solving  problem  PI,  and  have  demonstrated  that  it  is 
possible  to  significantly  reduce  the  computational  requirement  of  the 
Fox-Landi  method  by  solving  an  approximation  problem  to  obtain  a good 
estimate  of  an  appropriate  value  for  the  Lagrangian  multiplier.  In 
this  section  we  describe  a different  way  to  reduce  the  computational 
requirements  of  all  the  algorithms  that  have  been  discussed.  As  can 
be  seen  by  reexamining  Sherbrooke's  approach  (see  problem  P6)  and  the 
Fox-Landi  algorithm  (see  problems  P8  and  P9) , the  amount  of  computation 
required  to  solve  problem  PI  using  these  methods  is  directly  proportional 
to  the  number  of  depot  stock  levels  explicitly  examined.  Consequently, 
if  this  number  can  be  reduced,  then  the  total  time  required  to  compute 
an  optimal  solution  can  also  be  reduced. 

The  method  that  we  describe  in  this  section  to  estimate  the  optimal 
depot  stock  level  will  be  of  particular  value  when  the  expected  number 
of  units  in  the  depot  resupply  system  for  an  item  is  20  or  more.  The 
approximation  algorithm  can  reduce  computation  time  for  the  algorithms 
described  in  Sec.  IV  by  as  much  as  90  percent  for  these  high  demand 
items . 
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5.(s  ; sQ)  = l (x  - s )p (x | X T (sQ) ) , 

J J *S  p J J J 


for  s.  > 0 and  integer,  and  the  piecewise  linear  completion  of  B. , 

J A J 

call  it  B.,  as  follows: 

J 


IB j ( t ; Sq)  if  t is  a non-negative  integer. 

[Bj(Sj;  sQ)  - (sj  - 1;  sQ)](t  _ (sj  “ D) 

+ B(Sj  - 1;  sQ) , Sj  - 1 < t < Sj  ; 

where  s^  is  a non-negative  integer, 
and  B (—1 ; sQ)  = 00 . 


AB.(s.;  s0)  = B.(s.;  sQ)  - B.(s.  - 1;  sQ) 


when  s^  is  a non-negative  integer,  and 


D(s.;  s ) = {v:  AB . (s . ; sn)  < v < AB.(s.  + 1;  s„)}  . 

J 0 J J 0 1 j 0 


Observe  that  D.Cs^;  Sq)  U {ABj(s.;  Sq)}  is  the  set  of  subgradients  of 
Bj  at  Sj . Then  an  alternative  way  of  verifying  that  s*  is  an  optimal 
base  stock  level  is  to  show  that  - 0c  E D(s*;  Sq) . 

Next  let 


A 111 

V s2 ’ sn;  s0)  = A (Bj(Sj!  S0}  + 0CSj)  ' 

J = 1 


By  dropping  both  the  integrality  and  non-negativity  restrictions  on 
Sq,  we  obtain  the  following  relaxation  of  problem  P8: 

min  / 0cs„  + min  {F(s,,  ...,  s ; sn) : s fixed}?  . (P10) 

sn  1 ° s =0, 1, . . . 1 n 0 0 I 
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If  Sq  is  the  optimal  solution  to  problem  P10,  then 


3F 

3s„ 


+ 0c  = 0 . 


(E3) 


But 


3L-  v l 

3s0  ’ j-1  3Tj  Ss0 


Furthermore,  by  writing  B^  (s^  ; sq)  as 


l Kp(K  + s | XT  (s  ))  , 
K=1  3 J 


we  see  that 


K+s.-l 

SB.  « -X.T.(sn)  (X.T. (s„))  3 

-1  = l X . Ke  3 3 ° — LJ_Q_ 


3TJ  K=1  3 


(K  + s - 1) ! 


K+s . 

X4T,(sn)  (X.T.  (s„) ) 3 


l KX.e  3 3 0 


K=1 


(K  + s ) ! 


-Xj AB(Sj ; sQ)  . 


As  we  discussed  in  Sec.  II,  the  function 


B0(s0)  = l (x  ~ S0>P (x I XD) 


x>s 


can  be  closely  approximated  by  an  exponential  function  of  the  form 

_boso 

a^e  , where  a^  and  b^  are  positive  real  numbers.  Then 


-31- 


and 


VV  - vv  (1  - v V X v Vo 


!Li  (1  - V „ -boso 

X~  a0b0e 


Upon  combining  these  observations  we  see  that 


IU 

fiT~  V£(sr  so) 

0 j=l  J J 


(1  - ri}  . "boso 

x aoboe 


Recall  that  -0c  e D(s*;  s^).  Consequently  -0c  approximates  the 
marginal  reduction  in  backorders  at  base  j when  the  stock  level  at  that 
base  is  s*.  After  making  this  substitution  and  representing  this  fur- 
ther  approximation  of  9F/3sq  by  9F/8sq,  we  see  that 


I?  1 

~ — = - £ (1  - r.)X.  y 0ca_b  i 

9s0  j=l  J J x 00 


"b0S0 


= -9ca0b0e 


■boso 


Substituting  this  approximation  into  Eq.  (E3)  we  obtain  the  follow- 
ing estimate  of  the  optimal  depot  stock  level: 


s0  = - \ ~T~ 

bo  ( aobo 


(E4) 


Recall  the  value  of  §q  is  derived  based  on  an  exponential  approxi- 
mation of  ESq(Sq).  As  the  average  number  of  units  in  the  depot  repair 
cycle  increases,  that  is,  as  XO  increases,  the  quality  of  this  expo- 
nential approximation  improves  in  the  region  in  which  the  optimal  depot 
stock  level  should  be  located.  Consequently,  the  approximation  should 
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be  most  accurate  in  these  cases.  But,  the  problems  for  which  the 
search  for  the  optimal  depot  stock  level  is  most  time-consuming  for 
the  algorithms  described  in  Sec.  IV  correspond  to  the  items  having  a 
large  number  of  units  in  depot  repair.  Therefore,  the  proposed  ap- 
proximation method  will  be  most  appropriate  for  the  items  requiring 
the  greatest  amount  of  computational  effort. 

The  approach  we  have  described  for  estimating  the  optimal  depot 
stock  level  has  been  coded  and  tested  using  a sample  of  40  F-15  air- 
craft items.  The  test  consisted  of  two  sets  of  runs.  In  the  first  set, 
monthly  flying  was  divided  among  3 bases;  in  the  second  set  the  same 
monthly  flying  program  was  divided  among  5 bases.  The  total  budget 
distributed  among  the  40  items  ranged  from  $34  million  to  $65  million 
in  the  first  set  of  runs,  and  from  $34  million  to  $88  million  in  the 
second  set.  Table  3 contains  the  data  indicating  both  the  optimal  and 
estimated  depot  stock  level  for  each  item  in  both  runs. 

As  shown  in  the  table,  there  is  usually  no  single  optimal  depot 
stock  level  for  an  item.  Rather  the  optimal  value  depends  on  the 
amount  of  total  item  system  stock  available  for  distribution  among 
the  depot  and  bases.  The  estimate  of  optimal  depot  stock  is  quite 
close  to  the  optimal  value  in  all  cases.  Furthermore,  the  increase 
in  expected  system  backorders  using  the  estimated  depot  stock  levels 
rather  than  the  optimal  levels  is  generally  small.  For  most  items 
the  increase  is  substantially  less  than  .1  backorders. 

The  results  of  the  tests  indicate  that  it  is  possible  to  estimate 
closely  the  optimal  depot  stock  level  using  Eq . (E4).  Additionally, 
incorporating  this  method  for  estimating  the  optimal  depot  stock  into 
the  algorithms  described  in  Sec.  IV  will  considerably  reduce  the  search 
required  to  find  the  optimal  depot  stock  level  and  will  therefore 
markedly  reduce  the  computational  time  needed  to  solve  problem  PI 
using  these  algorithms. 
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